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ABSTRACT 

Context. Ices are present in comets and in the mantles of interstellar grains. Their chemical composition has been indirectly derived by 
observing molecules released in the gas phase, when comets approach the sun and when ice mantles are sublimated or destroyed, e.g. 
in the hot cores present in high-mass, star-forming regions. Comparison of these chemical compositions sheds light on the formation 
of comets and on the evolution of interstellar matter from the molecular cloud to a protoplanetary disk, and it shows, to first order, a 
good agreement between the cometary and interstellar abundances. However, a complex O-bearing organic molecule, ethylene glycol 
(CHjOHE, seems to depart front this correlation because it was not easily detected in the interstellar medium (Sgr B2) although it 
proved to be rather abundant with respect to other O-bearing species in comet C/1995 Ol (Hale-Bopp). Ethylene glycol thus appears, 
together with the closely related molecules glycolaldehyde CHiOHCHO and ethanol CH 3 CH 2 OH, as a key species in the comparison 
of interstellar and cometary ices as well as in any discussion on the formation of cometary matter. 

Aims. It is important to measure the molecular abundances in various hot cores to see if the observed differences between the interstel¬ 
lar medium and the comets are general. We focus here on the analysis of ethylene glycol in the nearest and best studied hot core-like 
region, Orion-KL. 

Methods. We use ALMA interferometric data because high spatial resolution observations allow us to reduce the line confusion 
problem with respect to single-dish observations since different molecules are expected to exhibit different spatial distributions. 
Furthermore, a large spectral bandwidth is needed because many individual transitions are required to securely detect large organic 
molecules. Confusion and continuum subtraction are major issues and have been handled with care. 

Results. We have detected the aGg’ conformer of ethylene glycol in Orion-KL. The emission is compact and peaks towards the Hot 
Core close to the main continuum peak, about 2" to the south-west; this distribution is notably different from other O-bearing species. 
Assuming optically thin lines and local thermodynamic equilibrium, we derive a rotational temperature of 145±30 K and a column 
density of 4.6±0.8 10 15 cnT 2 . The limit on the column density of the gGg’ conformer is five times lower. 

Key words. Astrochemistry - ISM: molecules - Radio lines: ISM - ISM: individual objects: Orion-KL - Comets: general 


1. Introduction 

The good global correlation between the chemical abundances 
in cometary ices and molecular hot cores (see Fig.3 of Bockelee- 
Morvan et al. 2000) is a strong argument in favour of similar for¬ 
mation processes in these two rather different types of objects. 
Ion-molecule and surface reactions at a low (< 200 K) tempera¬ 
ture are thus invoked (Herbst & van Dishoeck 2009) but one may 
not yet exclude the possibility of a very direct connection imply¬ 
ing comets as being "frozen pieces of the interstellar medium" as 
discussed as early as 1975 by Shimizu (see also Despois 1992, 
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and references therein, and Mumma & Charnley 2011, for a re¬ 
cent review). 

However, at least two important problems arise: the D/H ra¬ 
tio (e.g. Ceccarelli et al. 2014), and the crystalline/amorphous 
ratio in silicates do not seem to match well (e.g. Crovisier et al. 
1997). Ethylene glycol may be a third example of a discrepancy. 
This opens the possibility for a wide range of models, ranging 
from pure interstellar matter to pure products of nebular chem¬ 
istry, to explain the formation of comets. Consequently, differ¬ 
ent models could explain the protosolar nebula and solar system 
formation (see e.g. Irvine et al. 2000). This in turn has conse¬ 
quences on the way the collapse proceeds in young stellar ob¬ 
jects, on physical and chemical evolution and timescales for the 
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collapse and pre-collapse phase and later mixing in the disk due 
to turbulence (Bockelee-Morvan et al. 2002). 

Ethylene glycol, HOCH 2 -CH 2 OH, is a dialcohol, a molecule 
chemically related to standard alcohol (ethanol, CH 3 -CH 2 OH). 
Ethylene glycol is also known as an antifreeze coolant for car en¬ 
gines. Ethylene glycol was first found in the interstellar medium 
towards the Galactic centre along the line of sight to Sgr B2 de¬ 
spite line blending problems (Hollis et al. 2002; Requena-Torres 
et al. 2008). In comets, it has been identified for the first time 
in the Comet C/1995 Ol (Hale-Bopp) spectra by Crovisier et al. 
(2004a). 

The detailed comparison of the abundances of ethylene gly¬ 
col in the interstellar medium and in comets reveals a very dif¬ 
ferent situation: 

- In Hale-Bopp, ethylene glycol is rather abundant (0.25 % 
by number with respect to H 2 O, or 1 % by mass), just behind the 
main CHO species, methanol CH 3 OH (2.5 % by number), and 
comparable to formaldehyde (0.1-1.1 % by number). Ethylene 
glycol has recently been detected in two other comets: C/2012 F 6 
(Lemmon) and C/2013 R1 (Lovejoy), with similar abundances 
relative to water (0.24 and 0.33 %, Biver et al. 2014). The simpler 
ethanol is not detected (< 0.1 % - 3cr limit), nor the companion 
molecule to ethylene glycol, glycolaldehyde (HOCH 2 -CHO) (< 
0.04 %) (Crovisier et al. 2004b). 

- In the interstellar medium, ethylene glycol is not easy to 
identify, and until now it has been found towards only a few 
sources: Sgr B2 (Hollis et al. 2002; Requena-Torres et al. 2008), 
NGC 7129 FIRS 2 (Fuente et al. 2014), and NGC 1333-IRAS2A 
(Maury et al. 2014). Along the Sgr B2 line of sight, glycolalde¬ 
hyde is seen with a similar abundance (within a factor of 2, Hollis 
et al. 2002). Ethanol is more abundant by an order of magnitude 
in Sgr B2 (Nummelin et al. 2000), and seen in many sources. 

How can these differences be explained? Is the line of sight 
towards the galactic centre atypical, or are we seeing general 
differences between interstellar and cometary ices? Observations 
of hot cores, which appear to be the interstellar regions most 
similar to comets, are required to answer these questions. 

Orion-KL is one of the main targets to search for new 
molecules in the interstellar medium, because of its molecular 
wealth and its proximity. Interferometric maps are absolutely 
needed because they notably reduce the confusion problem due 
to line blends in the single dish spectra. As shown by Wright 
et al. (1996) and confirmed by our data (Guelin et al. 2008; 
Favre et al. 2011; Peng et al. 2012, 2013; Brouillet et al. 2013), 
different sets of molecules peak towards different positions in 
Orion-KL, and possibly even do not overlap on high (1-2") spa¬ 
tial resolution maps. A second clear advantage of interferometer 
maps, is the information they provide: if the spatial distribution 
is identical for various lines, they most likely arise from the same 
species and not from a single interloper exhibiting the same fre¬ 
quencies as those of the molecular species under study. This is 
especially true for ethylene glycol, a molecule with a complex 
microwave spectrum. 

In this paper, we report the results of our search for the two 
most stable conformers of ethylene glycol towards Orion-KL. 
The spatial distribution of the ethylene glycol emission and the 
derived physical parameters are presented in Sect. 3. In Sect. 
4 we discuss the chemistry in comparison to other interstellar 
sources and comets. We present the continuum emission map 
and its variation at nearby frequencies in Sect. 5. Finally our 
conclusions are summarized in Sect. 6 . 


2. Observations 

We primarily use here the ALMA Science Verification (ALMA- 
SV) data 1 taken on 20 January 2012 towards the IRc2 region in 
Orion. The observations were carried out with 16 x 12 m anten¬ 
nas in the Band 6 frequency range from 213.715 to 246.627 GHz. 
20 slightly overlapping spectral windows were observed, each 
with a 1.875 GHz bandwidth and 3 840 channels resulting in a 
0.488 MHz channel spacing corresponding to a velocity sepa¬ 
ration of 0.64 km s -1 . The observations were centred on coor¬ 
dinates (072000 = 05 /, 35'”14?35, 6 72000 = -05°22 , 35'.'00) with a 
primary beam of ~27". 

We used the CASA software 2 for initial processing, and then 
we exported visibilities to the GILDAS package 3 for further 
analysis. The line maps were cleaned using the CLARK algo¬ 
rithm (Clark 1980). The synthesized beam is 1.90" x 1.40" with 
a PA of 170° and the brightness temperature to flux density con¬ 
version factor is 9 K for 1 Jy per beam. The continuum emission 
was subtracted in the maps by carefully selecting line-free chan¬ 
nels (see Sect. 5). 

We also analyzed several observational data sets obtained 
with the IRAM Plateau de Bure Interferometer (see Table 1 in 
Favre et al. 2011). These observations cover about 2 GHz from 
80.5 GHz to 226 GHz with a velocity separation ranging from 
0.42 km s 1 to 2.33 km s 1 and a spatial resolution ranging from 
1.79" x 0.79" to 7.63" x 5.35", depending on the frequency 
and the IRAM interferometer configuration. Four data sets corre¬ 
spond to frequencies included in the ALMA-S V data. In particu¬ 
lar, the IRAM 223.4—223.9 GHz observations have spectral and 
spatial resolution and sensitivity similar to those in the ALMA 
data set for the same spectral range. 

3. Results 

3.1. Conformers of ethylene glycol and selection of lines 

Ethylene glycol is an asymmetric top molecule with coupled ro¬ 
tation around its C-C bond and its two C-O bonds leading to 
different spatial arrangements or conformers (see Christen et al. 
1995, 2001). There are six conformers relative to the C-C bond 
with gauche (G) arrangement of the two ethylene groups and 
four less stable conformers with anti-arrangement of the CH 2 
groups. In the G arrangement, the two OH groups adopt the 
gauche orientation with respect to each other. The notations a 
and g are used for the anti (a) and gauche (g) orientations of the 
OH groups with respect to the two C-O bonds. When the OH 
group provides the hydrogen intramolecular bonding, the gauche 
orientation is noted g’. There are only two conformers with in¬ 
tramolecular bonding, aGg’ and gGg’, and Muller & Christen 
(2004) showed that the latter lies about 200 cm -1 or about 290 K 
higher in energy than aGg’. For these two conformers, tunneling 
is observed between two equivalent equilibrium configurations 
and splits each rotational level into two distinct states designated 
v = 0 and v = 1 (e.g. Christen et al. 2001). The v = 1 state is 
around 7 and 1.4 GHz higher than the v = 0 state for the aGg’ 
and gGg’ conformers, respectively (Christen & Muller 2003). 

The rotation-tunnelling microwave spectra of the aGg’ and 
gGg ’ conformers were accurately determined in the 54-370 GHz 
(Christen & Muller 2003) and 77-579 GHz (Muller & Chris¬ 
ten 2004) ranges, respectively. The results from these works and 
earlier works are now available from the Cologne database for 

1 http://almascience.eso.org/almadata/sciver/OrionKLBand6/ 

2 http://casa.nrao.edu 

3 http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. Superposition of the unblended and partially blended ethylene 
glycol lines listed in Table 1. The lines corresponding to close pairs of 
transitions separated by 0.5-1 MHz are not included in this plot as their 
profile is broadened. The spectra are towards the ethylene glycol peak 
(a /2000 = 05 / '35 m 14547, Sj 2 ooo = -05°22'33'.'17) and are normalized to 
1 . 

molecular spectroscopy (CDMS 4 , Muller et al. 2001, 2005) and 
can be used for reliable identification of the ethylene glycol tran¬ 
sitions. Earlier ethylene glycol line predictions were used to as¬ 
sign to the aGg’ conformer five features observed by Hollis et al. 
(2002) in Sgr B2; two other features were also assigned to the 
same conformer in the same source by Requena-Torres et al. 
(2008). The CDMS was used to identify more than 10 aGg’ con¬ 
former lines in the Hale-Bopp comet (Crovisier et al. 2004a). 
We note that several ethylene glycol transitions observed in the 
direction of Sgr B2 are blended, while transitions are well ’iso¬ 
lated’ towards Hale-Bopp. 

We used the most recent CDMS data to search for transitions 
for the aGg’ and gGg’ conformers in the ALMA-SV data cube. 
From line intensity predictions at 120 K, a typical temperature 
for the region (e.g. Favre et al. 2011), we have examined the 
strongest lines among the ethylene glycol transitions expected in 
the different spectral bands. We produced synthetic spectra for 
both conformers (see Sect. 3.4 and 3.5). We also derived syn¬ 
thetic spectra for higher temperatures of the gGg ’ conformer to 
investigate the presence of the gGg ’ conformer in hotter gas be¬ 
cause gGg’ is higher in energy than aGg’ (see above). 

3.2. Spatial distribution and line profile of the aGg’ conformer 
of ethylene glycol 

In the AFMA dataset, most ethylene glycol lines suffer mild 
to heavy blending. Table 1 lists 57 transitions that are mildly 
blended and a noncontaminated transition lying at 228 752 MHz. 
All these transitions present a peak in their spatial distributions 
at («/2000 = 05 A 35 m 14?47, <5,2000 = -05°22'33'.'17). In Fig. 1 we 
plotted a superposition of the spectra of the partially blended 
ethylene glycol lines towards this position, excluding line dou¬ 
blets with a difference in frequency of ~ 0.5-1 MHz. We nor¬ 
malized the spectra to a peak intensity of 1. All line profiles are 
similar in the 5 to 10 km s 1 velocity range. They can be fitted 
with a Gaussian centred at a velocity of 7.6 km s“ 1 and a line 
width of 2.3 km s -1 . We have averaged the transitions that ex¬ 
hibit no blends around 5 km s 1 on one hand, and the transitions 
that exhibit no blends around 10 km s 1 on the other hand. Fig¬ 
ure 2 shows the superposition of these two spectra together with 
the 228 752 MHz transition spectrum. The resulting line profile 
of the unblended parts fits perfectly the profile of the 228 752 
MHz transition. The Gaussian fit shows a faint excess around 

4 http://www.astro.uni-koeln.de/cdms 



Fig. 2. Superposition of the 228 752 MHz Uansition spectrum at the 
ethylene glycol peak (in black) with the average spectrum of 11 transi¬ 
tions that have no blends near 5 km s -1 (in red) and the average spectrum 
of 14 transitions that have no blends near 10 km s -1 (in blue). All spec¬ 
tra are normalized to 1. A Gaussian fit to the 228 752 MHz line is shown 
in green. A faint wing is visible at lower velocities. 

5 km s -1 . This excess is significant as it is above the noise of 
the 228 752 MHz individual spectrum. Furthermore it appears 
on both profiles, which use different transitions, shown in Fig. 2. 

The spatial distribution of ethylene glycol in Orion-KF 
is shown in Fig. 3, which presents the intensity maps in¬ 
tegrated over the noncontaminated channels of the five best 
lines (transitions at 219 385 MHz, 222 348 MHz, 228 752 MHz, 
231524 MHz, and 235 834 MHz). The 228 752 MHz transition 
appears to be the only one free of contamination, i.e. without any 
blending line above the noise level; the channel maps are shown 
in Fig. 4. 

The aGg ’ conformer distribution is compact and appears to 
be slightly bigger than the beam. We find an upper limit for the 
deconvolved size of 2.4"x 1.1". The emission peak is located 
~ 2" to the south-west of the peak of continuum emission (see 
map in lower right panel of Fig. 3 and Fig. 8 map); it lies between 
the methyl formate MF2 and MF6 peaks of Favre et al. (2011). 
The line profile is similar to the methyl formate and ethanol line 
profiles at the same position, centred at a velocity of ~7.6 km s 1 
and with a faint wing around 5 km s _I . A detailed comparison of 
these O-bearing molecules will be performed in a forthcoming 
paper. 

3.3. Spectral confusion 

The spectral line density in Orion is high (e.g. 8.6 lines per 
100 MHz in the IRAM 30m survey of Tercero et al. 2010), com¬ 
parable to that in Sgr B2 (6-10 lines per 100 MHz, Friedel et al. 
2004; Halfen et al. 2006) for high-sensitivity observations. 

Snyder et al. (2005) have given the essential criteria for es¬ 
tablishing the identification of a new interstellar molecule. Even 
if ethylene glycol has already been detected in Sgr B2 and two 
other interstellar sources, its identification is not obvious because 
of the weakness of the lines and spectral confusion. The different 
criteria are: 

- The rest frequencies must be known with a high degree of 
accuracy from laboratory measurements. 

- The central velocity and line width must agree for all de¬ 
tected transitions. 

- Blended lines must be separated by at least half the line 
width (FWHM) of the weaker line, to be considered as resolved. 

- Beam dilution must be taken into account when comparing 
observations made at different frequencies. 
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Table 1 . Ethylene glycol transitions (aGg' conformer) not blended or only partially blended in the ALMA dataset. 


Frequency a 
(MHz) 

Quantum numbers b 


Sp 2c 

(D 2 ) 

log(A ul ) " 

E d 
(K) 

£ d 
Oil 

214526.41410 

19(11, 8) v = 1-19(10, 9) v = 1 

50.7 

-4.67073 

153.1 

273 

214526.41470 

19(11, 9) v = 1 - 19(10,10) v = 

1 

65.2 

-4.67067 

153.1 

351 

214712.96650 

21(5,16) v = 0-20(5,15) v= 1 


732.5 

-3.66141 

127.1 

387 

214 808.01940 

23(2,22) v = 0-22(2,21) v= 1 


685.7 

-3.61897 

132.4 

329 

214 829.92520 

23(1,22) v = 0 - 22(1,21) v= 1 


875.5 

-3.62187 

132.4 

423 

216 685.81470* 

21(3,19) v= 1 - 20(3,18) v = 0 


663.4 

-3.69253 

117.3 

387 

216 826.11240 

20(5,15) v= 1 - 19(5,14) v = 0 


557.3 

-3.74666 

116.8 

369 

217139.72350 

21(4,17) v = 0-20(4,16) v= 1 


786.8 

-3.61568 

123.9 

387 

217449.99470 

24(1,24) v= 1 - 23(1,23) v= 1 


928.6 

-3.59857 

136.4 

441 

217450.27020* 

24(0,24) v = 0 - 23(0,23) v = 1 


722.2 

-3.59863 

136.4 

343 

217 587.54780* 

21(2,19) v= 1-20(2,18) v = 0 


666.2 

-3.57614 

117.2 

301 

219 089.72000 

22(10,13) v = 0 - 21(10,12) ) v : 

= 1 

670.7 

-3.69311 

173.5 

405 

219 089.72790 

22 (10,12) v = 0 - 21(10,11) v = 

1 

521.6 

-3.69317 

173.5 

315 

219 384.91030* 

26(11,16) v = 1- 26(10,16) v = 

0 

15.8 

-5.28218 

232.4 

371 

219385.17780* 

22(9,14) v = 0-21( 9,13) v= 1 

703.8 

-3.67045 

164.3 

405 

219 385.32400* 

26(11,15) v = 1-26(10,17) v = 

0 

20.3 

-5.28213 

232.4 

477 

219 385.42560* 

22( 9,13) v = 0 - 21( 9,12) v = 1 

547.3 

-3.67051 

164.3 

315 

221007.82310 

21(4,18) v= 1-20(4,17) v = 0 


818.3 

-3.57562 

122.1 

387 

221038.79970* 

22(6,17) v = 0-21(6,16) v= 1 


771.6 

-3.62074 

142.6 

405 

221 100.31520 

22(5,18) v = 0-21(5,17) v= 1 


762.5 

-3.62553 

137.4 

405 

222 054.23240* 

21 (10,12) v= 1 -20(10,11) v = 

0 

624.5 

-3.68686 

163.0 

387 

222 054.23580* 

21 (10,11) v = 1-20(10,10) v = 

0 

485.8 

-3.68681 

163.0 

301 

222 348.60040* 

23(1,23) v= 1 - 22(1,22) v = 0 


889.5 

-3.57018 

126.0 

423 

222 349.15090* 

23(0,23) v = 1 - 22(0,22) v = 0 


692.0 

-3.57003 

126.0 

329 

223 741.66380* 

21(6,16) v= 1-20(6,15) v = 0 


738.2 

-3.60437 

132.0 

387 

226 643.30070* 

25(1,25) v = 0-24(1,24) v= 1 


752.5 

-3.54420 

147.7 

357 

226 643.45630* 

25(0,25) v = 0 - 24(0,24) v = 1 


967.6 

-3.54415 

147.7 

459 

228 602.83320 

23(16,7) v = 0-22(16,6) v= 1 


456.4 

-3.82381 

261.3 

423 

228 602.83320 

23(16,8) v = 0-22(16,7) v= 1 


354.9 

-3.82387 

261.3 

329 

228 752.77650*' 

23(2,22) v= 1 - 22(2,21) v = 0 


875.2 

-3.54020 

132.7 

423 

229 233.71610 

23(11,13) v = 0-22(11,12) v = 

1 

530.2 

-3.64599 

195.1 

329 

229 233.71660 

23(11,12) v = 0-22(11,11) v = 

1 

681.8 

-3.64594 

195.1 

423 

231 127.40080* 

23(7,16) v = 0-22(7,15) v= 1 


801.1 

-3.56515 

160.2 

423 

231524.03310* 

23(6,18) v = 0-22( 6,17) v= 1 


618.1 

-3.56638 

154.1 

329 

231564.00490 

24(1,24) v= 1 - 23(1,23) v = 0 


722.1 

-3.51672 

136.8 

343 

231564.31950 

24(0,24) v = 1 - 23(0,23) v = 0 


928.6 

-3.51666 

136.8 

441 

232 095.73820 

22 (12,11) v= 1-21(12,10) v = 

0 

462.3 

-3.67043 

195.4 

315 

232 095.73830 

22(12,10) v= 1 - 21(12,9) v = 0 

594.5 

-3.67038 

195.4 

405 

232 349.45090 

22 (20,2) v= 1 -21(20,1) v = 0 


146.9 

-4.27619 

320.9 

405 

232 349.45090 

22(20,3) v = 1 - 21(20,2) v = 0 


114.2 

-4.27614 

320.9 

315 

232 350.05940 

22(10,13) v= 1 -21(10,12) v = 

0 

522.1 

-3.61615 

173.8 

315 

232 350.06810 

22 (10,12) v = 1 -21(10,11) v = 

0 

671.4 

-3.61610 

173.8 

405 

232 881.53340 

23(6,17) v = 0-22(6,16) v= 1 


750.3 

-3.58376 

154.2 

423 

235 304.05030* 

22(6,16) v= 1 - 21(6,15) v = 0 


773.7 

-3.53805 

143.1 

405 

235 441.43720 

10(7,4) v = 0 - 9(6,3) v = 0 


56.6 

-4.34165 

51.1 

189 

235 441.48820 

10(7,3) v = 0 - 9(6,4) v = 0 


44.1 

-4.34161 

51.1 

147 

235 442.34320 

10(7,4) v = 1 - 9(6,3) v = 1 


42.1 

-4.36183 

51.4 

147 

235 442.39560 

10(7,3) v = 1 - 9(6,4) v = 1 


54.1 

-4.36177 

51.4 

189 

235 471.89630 

8(8,0) v = 0 - 7(7,1) v = 0 


49.0 

-4.20372 

49.1 

119 

235 471.89630 

8(8,1) v = 0 - 7(7,0) v = 0 


63.0 

-4.20376 

49.1 

153 

235 834.23970* 

26(1,26) v = 0-25(1,25) v= 1 


1 006.5 

-3.49194 

159.3 

477 

235 834.32720* 

26(0,26) v = 0 - 25(0,25) v = 1 


782.7 

-3.49199 

159.3 

371 

239 605.33240* 

24(11,14) v = 0-23(11,13) v = 

1 

728.5 

-3.57758 

206.9 

441 

239 605.33370* 

24(11,13) v = 0-23(11,12) v = 

1 

566.7 

-3.57753 

206.9 

343 

241545.26260 

24(7,18) v = 0-23(7,17) v= 1 


841.5 

-3.50446 

172.1 

441 

242 947.99050* 

23(9,15) v= 1 - 22(9,14) v = 0 


749.2 

-3.52928 

175.9 

423 

242 948.59120* 

23(9,14) v= 1 - 22(9,13) v = 0 


582.8 

-3.52923 

175.9 

329 

243 259.73980 

24(6,18) v = 0-23(6,17)v = 0 


324.4 

-3.80007 

166.2 

343 

244685.14770 

24(2,22) v= 1-23(2,21) v= 1 


939.2 

-3.43992 

150.1 

441 


a Transitions used for the rotational diagram (see Fig. 5) are indicated by a star. 

b The quantum numbers are J(K. d ,K c ) as defined in the CDMS database where v is a state number (see Sect. 3.1). 
c S is the line strength and p the dipole moment. 
d See Sect. 3.4 for the definition of the parameters. 
e Not blended transition. 
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Fig. 3. ALMA ethylene glycol integrated intensity maps of the five less 
blended lines. The frequency of the transition is indicated at the bot¬ 
tom of the map and the first contour level corresponds to 4<x. The beam 
shown in the bottom left corner is 1.9" x 1.39". The BN object posi¬ 
tion is (or72000 = 05 ,1 35"'14il 1, d/2000 = -05°22'22'. / 72), the radio source 
I position is (a J2 ooo = 05*35' B 14!51, d/2000 = -05°22'30'.'57), and the 
IR source n position is (a/2000 = 05 ,1 35"’14?36, d/2000 = -05°22'32'.'72) 
(Goddi et al. 2011b). The map to the lower right is a blow up of the 
228 752 MHz map. The triangle marks the position of the continuum 
emission peak (see map in Fig. 8) and the stars mark the MF2 and 
MF6 methyl formate positions of Favre et al. (2011). The ethylene gly¬ 
col emission peaks at (a/2000 = 05 A 35 m 14547, d/2000 = -05°22'33717). 
The lines are integrated over the velocity range 5-10 km s _1 for the 
222 348 MHz transition, 5.5-10 km s' 1 for 235 834 MHz, 7-10 km s' 1 
for 231 524 MHz, 6-9.5 km s -1 for 219 385 MHz and 4-11 km s' 1 for 
228 752 MHz. The level step and first contour are 6 and 2.8 K km s' 1 for 
the 222 348 MHz transition, 7.6 and 3.6 K km s -1 for 235 834 MHz, 2.5 
and 1.7 K km s” 1 for 231 524 MHz, 5 and 1.9 K km s“‘ for 219 385 MHz 
and 4.4 and 2.9 K km s" 1 for 228 752 MHz. 


- The relative intensities of the different transitions must be 
consistent (when, as we suppose it is the case here, an excitation 
close to local thermodynamic equilibrium is likely). 

- All transitions with intensity predictions leading to de¬ 
tectable signal levels must be present. 

For interferometric observations, another criterion can be 
added: the coherence of the spatial distribution. Indeed the spa¬ 
tial distribution of the different transitions must be similar for 
similar excitation conditions (e.g. for a similar upper level en¬ 
ergy). The probability of chance alignment in frequency in one 
spectrum may be high, but a chance alignment for all spectra in 
an extended map is unlikely. As most of the molecules display 
different distributions on a high-resolution scale, confusion will 
not produce similar maps for every transition. Requiring a sin¬ 
gle interloper for different transitions makes a misidentification 
much less probable. 
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Fig. 4. ALMA channel maps of the ethylene glycol 228 752 MHz transi¬ 
tion which has no apparent contamination. The level step and first con¬ 
tour are 0.44 K. The LSR velocity in km s _1 is indicated in the upper 
left corner of each map. 


Detailed study of the ethylene glycol spectroscopy is rela¬ 
tively recent and we can be confident that the rest frequencies 
are accurate (see Sect. 3.1). In our study of Orion-KL , there is a 
frequency agreement among all the partially blended transitions: 
the line centre velocities are in good agreement within the un¬ 
certainties (7.6+0.2 km s 1 ) as well as the line widths and the 
general profiles (see Fig. 1). In addition, the noncontaminated 
parts of these lines display the same spatial distribution (see Fig. 
3). The lines that we consider partially blended with others are 
separated by at least their half-maximum intensity. The relative 
intensities of the transitions are consistent: no expected strong 
lines are missing and a rotational temperature can be estimated 
based on these transitions (see Sect. 3.4). 

Halfen et al. (2006) add that the evidence for emission from 
a given molecule must be over a substantial frequency range and 
that an extensive, self-consistent data set is necessary to identify 
large organic species in the interstellar gas. Indeed, our initial 
search for ethylene glycol in Orion-KL based on limited spectral 
ranges with the IRAM Plateau de Bure Interferometer (PdBI), 
was not convincing enough. The ALMA Science Verification 
data are similar to our 223 GHz PdBI data for their sensitivity 
and their spectral and spatial resolution but they cover 32.9 GHz 
instead of 0.5 GHz. Using all the PdBI data sets (Favre et al. 
2011), which cover ~2 GHz, we have identified seven blended 
lines coinciding with aGg’ conformer transitions 5 . However, we 
have discarded these lines from the present analysis because of 
line blending problems and because of different spatial and spec¬ 
tral resolution used in our PdBI data. 


3.4. Physical parameters of the aGg’ conformer of ethylene 
glycol 

Assuming optically thin lines and local thermodynamic equilib¬ 
rium, we have estimated the rotational temperature and the col¬ 
umn density of ethylene glycol using the equation (e.g. Gold- 


? Frequencies of the PdBI ethylene glycol blended lines: 
105.7010 GHz, 105.8347 GHz, 203.4309 GHz. 203.4533 GHz, 
223.4589 GHz, 223.5205 GHz, 223.7417 GHz. 
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Fig. 5. Rotational diagram at the emission peak based on the less- 
contaminated ethylene glycol lines indicated by a star in Table 1. 


smith & Langer 1999): 


ln( 


8nk B v 2 W 

T1 —) = ln ^ 
hc J g u A al 


N 

QTT) 


)- 



(1) 


where W is the integrated line intensity (K km s 1 ), v the line fre¬ 
quency (Hz), A u i the Einstein coefficient for spontaneous emis¬ 
sion (s' 1 ), g u the degeneracy of the upper state, N the total col¬ 
umn density, Q(T) the partition function, E u the upper state en¬ 
ergy, and T the excitation temperature. 

The rotational diagram is obtained at the emission peak from 
the less-contaminated ethylene glycol lines to derive reliable 
physical parameters (see Fig. 5). Taking the partition function 
Q(T) = 234500 from the nonlinear plot of the CDMS values 
of Q at 145 K, we obtain a rotational temperature of 145+30 K 
and derive a total column density averaged over the synthesized 
beam N a Gg' of 4.6+0.8 10 15 cm' 2 . From the continuum map (see 
Sect. 5), we derived the beam averaged projected density N H , as 
in Favre et al. (2011), and we find an ethylene glycol abundance 
of ~10' 9 . 

We calculated the synthetic spectrum for a velocity v+sr of 
7.6 km s' 1 , a linewidth Avi /2 of 2.3 km s' 1 , a column density 
N a cg’ °f 4.6 10 15 cm' 2 and a temperature T of 145 K. Figure 6 
shows the synthetic spectrum overlaid on the observed ALMA 
spectrum towards the ethylene glycol peak in the range 228.74— 
228.79 GHz. Figure 7 shows the observed and synthetic spectra 
of the ethylene glycol transitions listed in Table 1. The synthetic 
spectrum overlaid to the whole ALMA band spectrum is shown 
in Fig. A.l. Figure 6 shows that the transition at 228 752 MHz 
is not blended, while the two other lines around 228 778 and 
228 782 MHz show line contamination and are not included in 
Table 1. 


3.5. Search for the gGg’ conformer of ethylene glycol 

Because the gGg’ conformer of ethylene glycol is not yet de¬ 
tected in the interstellar medium, we also searched for transi¬ 
tions from this stable conformer although it is 200 cm' 1 (290 K) 
higher in energy than the aGg’ conformer (Muller & Christen 
2004). 

All expected strong lines appear to be blended in the ALMA 
data cube used in this study. Making the simplifying, but admit¬ 
tedly questionable, hypotheses that the gGg ’ conformer presents 
the same spatial distribution and the same temperature as the 
aGg’ conformer, the comparison of the gGg’ synthetic spectrum 
with the ALMA observed spectrum at the aGg ’ peak leads to a 



228740 228750 228760 228770 228780 228790 

Rest Frequency (MHz) 

Fig. 6 . Part of the ALMA spectrum (in black) towards the ethylene 
glycol peak. The synthetic spectrum of ethylene glycol (EG) is over¬ 
laid in red. The parameters used for the synthetic spectrum are: Vlsr = 
7.6 km s' 1 , Avi /2 = 2.3 km s' 1 , N a a g ' = 4.6 10 15 cm' 2 , T = 145 K. Only 
the line on the left is devoid of confusion. 
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Fig. 7. ALMA spectrum (in black) towards the ethylene glycol peak at 
the ethylene glycol frequencies listed in Table 1. The synthetic spectrum 
of ethylene glycol (EG) is overlaid in red and the value of the upper state 
energy (E„) is indicated in blue. 


maximum column density of 10 15 cm' 2 , i.e. <1/5 of the aGg’ 
conformer. We have further assumed here that the partition func¬ 
tions of both conformers are identical since their rotational con¬ 
stants are very similar (Muller & Christen 2004). 

Because the gGg’ conformer is 290 K higher in energy than 
aGg ’ and thus may probe hotter gas, we have also derived a syn¬ 
thetic spectrum at the aGg’ peak for a temperature of 250 K 
significantly above the 145 K used in the previous Section. We 
find a projected density upper limit of 2.3 10 15 cm' 2 . With this 
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higher temperature, we have identified a few lines correspond¬ 
ing to gGg’ transitions 6 and with no apparent strong blending. 
However, the lines are too few and weak to dispel the confusion 
problem. 

The tentative detection of ethylene glycol towards 
IRAS 16293-2422, source B, made by Jprgensen et al. 
(2012) corresponds to the 220249.8 MHz transition of the 
gGg’ conformer. The narrow lines detected towards the IRAS 
solar-type source B by Jorgensen et al. show little or no line 
confusion. The presence of ethylene glycol, however, must still 
await confirmation from detection of other transitions especially 
those from the most stable aGg’ conformer. Unfortunately, 
the 220249.8 MHz gGg’ ttansition is partially blended in our 
ALMA spectrum towards Orion-KL . 


4. Discussion 

4. 1. Comparison with interstellar sources and comets 

First, we summarize our findings in Orion-KL. The ethylene gly¬ 
col peak is close to the prominent methyl formate position MF2 
of Favre et al. (2011) and the derived temperature for ethylene 
glycol is similar to that derived for methyl formate. However, 
the column densities obtained for ethylene glycol, around 4.6 X 

10 15 cm~ 2 , are not very sensitive to the rotational temperature in 
this temperature range (see Fig. B .6 of Peng et al. 2012, in the 
case of methanol). Favre et al. (2011) find an upper limit for the 
glycolaldehyde column density of 3.5 X 10 14 cm ~ 2 towards MF2 
for a temperature of 140 K and a spatial resolution of 3'.'8 X 270 
(line confusion does not allow us to get a significant upper limit 
at a higher spatial resolution). Peng et al. (2012) derive a column 
density of 4.2 + 0.4 X 10 18 cuT 2 for the methanol peak named 
dM-1 (same position as MF2) with the same conditions. Using 
the ALMA data, we have also studied the ethanol emission (De- 
spois et al. in prep.) and we find a column density of ~ 2.5 X 

10 16 cuT 2 and a temperature of ~ 140 K towards the ethanol 
peak of emission close to the ethylene glycol peak. 

Clearly, and as in Sgr B2 (Belloche et al. 2013) and 
NGC 7129 FIRS 2 (Fuente et al. 2014), ethanol is more abun¬ 
dant than ethylene glycol (by a factor of 5 in our data). But ethy¬ 
lene glycol is at least 10 times more abundant than glycolalde¬ 
hyde in Orion-KL whereas it is only twice more abundant in 
NGC 7129 FIRS 2 and both molecules have similar abundances 
towards Sgr B2. In all cases methanol is widespread and highly 
abundant. 

Ethylene glycol is also more abundant than glycolaldehyde 
in comets and less abundant than methanol (by a factor 7-10, 
Crovisier et al. 2004b; Biver et al. 2014), whereas ethanol has 
not yet been detected with an upper limit 2.5 times lower than 
ethylene glycol. 

Table 2 compares the column densities derived for ethy¬ 
lene glycol and related molecules towards Orion-KL, two other 
prominent interstellar sources and the 3 comets where ethylene 
glycol has been identified. We note that the comparison suffers 
from the fact that the spatial distribution is different for the differ¬ 
ent molecules and the spatial resolution is also different between 
the observed interstellar sources. 


6 Frequencies of gGg ’ conformer transitions corresponding to not very 
blended lines in the ALMA spectrum at the aGg’ peak: 220.2498 GHz, 
224.8281 GHz, 226.4111 GHz, 227.2706 GHz, 228.9154 GHz, 
228.9179 GHz, 228.9374 GHz, 228.9399 GHz, 229.4537 GHz, 
229.4540 GHz,230.4622 GHz, 232.7066 GHz, 232.7076 GHz, 
232.7080 GHz, 241.6567 GHz. 


With the recent detection of ethylene glycol in two comets 
(Biver et al. 2014), whereas more abundant interstellar species 
like ethanol and glycolaldehyde are still undetected, the problem 
of ethylene glycol formation gets even more acute. Pathways to 
its formation are still under investigation. Walsh et al. (2014) 
summarize available grain-surface scenarios which include a) 
formation from CH 2 OH radical, b) hydrogenation of OCCHO 
itself resulting from atomic oxygen addition on H 2 CO, c) some 
non detailed UV or CR processing of ices, d) possibility of ther¬ 
mal processing. 

4.2. Orion-KL 

Comparison of the ethylene glycol disttibution in Orion-KL with 
other available physical or chemical ttacers may shed some 
light on the molecular formation conditions in Orion-KL. We 
compared the maps of these ttacers using the Aladin software 7 
(Bonnarel et al. 2000). This in turn could help us to understand 
the molecular complexity in other regions and perhaps in comets. 

There is no obvious role played by temperature: the ethy¬ 
lene glycol peak does not show up in any of the high-resolution 
temperature maps of CH 3 OH (Beuther et al. 2006; Friedel & 
Widicus Weaver 2012) and NH 3 (Goddi et al. 2011a) or in any 
Mid-IR peaks (Okumura et al. 2011). 

The various ttacers do not show any concentration of mat¬ 
ter at the ethylene glycol peak. On the contrary, high-resolution 
maps of 13 CO (Zapata et al. 2011) and continuum emission 
present a hole at the ethylene glycol peak and no clear associ¬ 
ation is seen with NH 3 (Wilson et al. 2000). 

At higher spatial resolution two continuum spots, C21 and 
C22 (Friedel & Widicus Weaver 2011), appear on each side of 
the ethylene glycol peak. The northern spot is close to the SMA1 
continuum source and coincides with a 22 GHz H 2 O maser 
(Gaume et al. 1998) and the MF 6 methyl formate position of 
Favre et al. (2011). There is another 22 GHz maser spot in the 
south at the same position as the MF2 methyl formate peak near 
but not exactly coincident with the continuum peak. 

In the near-IR continuum as well as in H 2 2.2 pm (Lacombe 
et al. 2004, D. Rouan, private communication), the region ap¬ 
pears as a dark patch; this is also the case at longer wavelength 
(e.g. at 20 /rm, Robberto et al. 2005). 

One short CO jet among jets identified by Zapata et al. 
(2009) (and coming from the common "explosion centre" in be¬ 
tween source I and BN) overlaps the ethylene glycol peak. 

A precise comparison with ALMA maps of other oxygenated 
species is in preparation. According to previous studies by us 
and other groups, it appears clearly that the ethylene glycol dis¬ 
tribution is definitely different from that of methyl formate and 
dimethyl ether (see e.g. Brouillet et al. 2013). It is also differ¬ 
ent from CN bearing species, such as C 2 H 5 CN, from acetone 
CH 3 COCH 3 (Peng et al. 2013; Friedel & Snyder 2008; Widicus 
Weaver & Friedel 2012), and from deuterated methanol (Peng 
et al. 2012). One of the 13 C methanol peaks in Peng et al. (2012) 
is however coincident with the ethylene glycol peak. Compared 
with HDO, which is typically released from grain surfaces, ethy¬ 
lene glycol appears to come from the same region as the main 
HDO emission delineated in the maps of Neill et al. (2013) ob¬ 
tained from the ALMA-SV data. However, looking at the indi¬ 
vidual channel maps of the two molecules, we find at each ve¬ 
locity that the ethylene glycol emission peak is ~ 1 " west of the 
HDO peak. 

7 Centre de Donnees Astronomiques de Strasbourg, http://aladin.u- 
strasbg. fr/aladin. gml 
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These first comparisons of maps: a) seem to exclude that the 
temperature deduced from present-day observations could play 
a dominant role in the observed differences (but an earlier gas 
warming at ethylene glycol position is not excluded); and b) 
show that ethylene glycol is distinct from species like dimethyl 
ether and methyl formate, which are strong in the compact ridge 
and likely to be produced on grain mantles. The compact nature 
of the ethylene glycol emission compared to the more extended 
distributions of methanol and ethanol calls for a special produc¬ 
tion scenario. The above elements suggest that one possibility 
could be a collision of a high-velocity jet identified in CO (but 
possibly ionized) with a clump at the edge or in front of the main 
hot core condensation, which would be responsible for the ethy¬ 
lene glycol production and the hole in continuum emission. This 
collision could also be favourable to the excitation of the two 
22 GHz maser spots. Another scenario, involving UV photons 
and reactions with radicals, is not excluded. However, there is 
no special observational evidence in favour of this scenario, as 
no local source emitting strong UV emission has been identified 
yet in this area, and penetration of external UV photons into such 
a dense region would require voids in the dust distribution. The 
association of two CH 2 OH radicals to form (CH 20 H) 2 , although 
appealing by its simplicity, is not favoured as the CH 2 OH radical 
has a low mobility due to OH bonds (Walsh et al. 2014). 

5. Continuum emission 

A map of the continuum emission at 1.3 mm in Orion-KL is a by¬ 
product of our efforts to isolate the line emission by subtracting 
the underlying continuum. The ALMA data covers a total band¬ 
width of 32.9 GHz centred at 230.171 GHz. For practical rea¬ 
sons, we have separated each one of the 20 spectral windows into 
four 1 000 -channel sub-windows, and we have searched for line- 
free channels in each sub-window using the following method. 

If a channel map has no noticeable line emission, the dif¬ 
ference between the maximum and the minimum intensities of 
the channel map will be lower than for channels with stronger 
molecular emission. Sorting out channels by increasing value of 
this difference, we have retained the best channel candidates to 
be line-free. 

We then used the fact that if two channel maps are devoid 
(within the noise) of molecular emission, their difference is noise 
plus the difference in continuum emission. This contribution is 
negligible if the channel frequencies are close enough (as is the 
case in a sub-window). We selected the best of channel pairs with 
a low and spatially uniform noise difference. 

Next, we eliminated channel pairs whose frequencies are 
closer than a typical line width, and might thus have identical 
molecular emission, inappropriately passing our second crite¬ 
rion. 

Then we obtained the continuum map for a given sub¬ 
window by averaging the two channel maps of the best pair. 

Finally, we checked the coherence of the continuum maps 
obtained for adjacent sub-windows. 

It was not possible to identify line-free channels for all the 
sub-windows, in that case, we took the continuum emission 
from a sub-window adjacent in frequency. We have thus ob¬ 
tained an 80-channel continuum cube in the range 213.7 GHz 
to 246.6 GHz. 

Figure 8 compares the continuum maps obtained at the same 
frequency and with a similar spatial resolution with ALMA (this 
work) and the IRAM Plateau de Bure Interferometer (Favre et al. 
2011). There is a good agreement between both maps which 
show the same four main components. The ALMA Hot Core 
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Fig. 9. Flux density per beam (in Jy/beam) versus frequency across 
ALMA band 6. Left: at the position of the ethylene glycol peak. Right: 
at the Hot Core position. 

position is (0+2000 = 05 / '35'"14?55, d/2000 = -05°22'3173) and it 
corresponds, within 0.2", to the Ca continuum peak identified in 
Favre et al. (2011). The slight differences observed in the ALMA 
and IRAM interferometer maps are likely due to different spatial 
filtering of the two arrays. 

Figure 9 shows the continuum variation with respect to the 
frequency over the whole ALMA-SV bandwidth. The spectrum 
is flat at the position of the ethylene glycol peak whereas it rises 
with frequency at the Hot Core position. This seems to be related 
to the sizes of the different continuum sub-sources associated 
with our main continuum component. In particular, if the contin¬ 
uum sub-source close to the Hot Core position is more compact 
than the synthesized beam, its intensity will vary as the beam 
size varies across the frequency range covered by the ALMA- 
SV data. We further note that our unresolved main continuum 
component (Hot Core) shows a NE-SW extension which encom¬ 
passes the two sub-sources seen in the sub-arcsec 870 pm dust 
continuum map of Tang et al. (2010) and the complex emission 
recently mapped with 0.5" ALMA resolution by Hirota et al. 
(2015). 

6. Conclusions 

1. We detected the aGg’ conformer of ethylene glycol in 
Orion KL. The emission peaks towards the Hot Core close 
to the main continuum peak, about 2 " to the south-west. 
The distribution is compact with a deconvolved size < 2.4"x 
1.1". It is clearly different from that of other complex O- 
bearing species like methyl formate, dimethyl ether, acetone, 
and ethanol. 

2. We deduced a beam averaged column density N a c g • of 
4.6+0 .8 10 15 ctrT 2 and a rotational temperature T of 
145+30 K from a population diagram. Using the continuum 
map, we find an ethylene glycol abundance with respect to 
H 2 of the order of 10 9 . 

3. Taking methanol CH 3 OH as a reference, the relative abun¬ 
dance of ethylene glycol is of the order of 10 ~ 3 in Orion- 
KL, similar to that in Sgr B2. We confirm with Orion-KL 
that the ethylene glycol abundance in interstellar sources is 
much lower than in the three comets where it has been de¬ 
tected with a ratio in the range 0.1-0.15. This result places 
ethylene glycol clearly outside the correlation found by 
Bockelee-Morvan et al. (2000) between Comet C/1995 Ol 
(Hale-Bopp) and four well-studied compact sources, includ¬ 
ing Orion-KL. 

4. The simple alcohol ethanol CH 3 CH 2 OH is more abundant 
than the di-alcohol CH 2 OHCH 2 OH by a factor of 5 in 
Orion-KL and even larger factors were found in Sgr B2 and 
NGC 7129. This result confirms a marked difference with the 
comets where, on the opposite, an upper limit of 0.5 has been 
derived for the ethanol/ethylene glycol ratio. 

5. We only identified weak lines of the gGg’ conformer at 
the expected frequencies, and all strong lines were blended. 
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Fig. 8 . Left: continuum map obtained with ALMA at 223 GHz. The beam shown in the bottom left comer is 1.92" x 1.38". Right: continuum 
map obtained with the IRAM Plateau de Bure Interferometer at 223 GHz (Favre et al. 2011). The beam shown in the bottom left comer is 1.79" x 
0.79". The level step and first contour for the two maps are 1 K. 


Hence we conclude that the gGg ’ conformer is not detected 
in our study, and we derive a column density less than 0.2 
times the aGg’ conformer’s at the emission peak of the lat¬ 
ter. 

6. The identification of the ethylene glycol aGg ’ conformer was 
made possible thanks to the ALMA-SV data. A large fre¬ 
quency range is necessary to identify complex molecules and 
interferometric maps of the transitions considerably helped 
us in diminishing the line confusion problem. 

7. Molecules have different spatial distribution even when a pri¬ 
ori of a same group (e.g. large O-bearing molecules). This 
should be taken into account in future comparisons between 
interstellar sources observed at different spatial resolution 
as well as between interstellar sources and comets. It might 
bring diversity in forming planets and their content in prebi- 
otic species. 
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Table 2. Comparison of the column densities and abundances of ethylene glycol and related molecules in interstellar sources and 
comets. 


Molecule 

Column density" 

(cm -2 ) 

Orion-KL SgrB2 NGC 7129 

C/1995 Ol 
(Hale-Bopp) 

Abundance" 
(%H 2 0) 
C/2012 F 6 
(Lemmon) 

C/2013 A1 
(Lovejoy) 

(CH 2 OH ) 2 

4.6 10 15 

2.3 10 1S 

2.0 10 15 

0.25 

0.24 

0.35 

CH 3 OH 

4.2 10 18 

5.6 10 18 

3.4 10 20 

2.4 

1.6 

2.6 

CH 3 CH 2 OH 

2.5 10 16 

2.1 10 16 

3.0 10 18 

< 0.1 

< 0.11 

<0.14 

CH 2 OHCHO 

<3.5 10 14 

OO 

O 

Ul 

1.0 10 15 

<0.04 

<0.08 

<0.07 

Relative abundance with respect to methanol 

(CH 2 OH ) 2 

1.1 10~ 3 

0.4 10 ~ 3 

0.006 10- 3 

0.1 

0.15 

0.13 

CH 3 OH 

1 

1 

1 

1 

1 

1 

CH 3 CH 2 OH 

6 10- 3 

4 10 ^ 3 

9 10 - 3 

<0.04 

<0.07 

<0.05 

CH 2 OHCHO 

<0.08 10- 3 

0.3 10 - 3 

0.003 10 - 3 

< 0.02 

<0.05 

<0.03 


a See Sect. 4.1 for references. 
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Fig. A.l. ALMA spectrum (in black) towards the ethylene glycol peak. The synthetic spectrum is overlaid in red. The parameters used for the 
synthetic spectrum are: Vlsr = 7.6 km s -1 , Avi/i = 2.3 km s -1 , N a c g ’ = 4.6 10 15 cm -2 , T = 145 K. Each box corresponds to a 120 MHz bandwidth 
and the ordinates are the intensities in K. The frequencies (in MHz) increase from left to right and from bottom to top. The transitions listed in 
Table 1 are indicated on the plots by "EG". 
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